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We show that a gauge singlet scalar S, with a coupling to the Higgs dou- 
blet of the form XsS^ SH^ H and with the S mass entirely generated by the 



Qh' Higgs expectation value, has a thermally generated relic density ^Is ~ 0-3 if 

■ ms ^ (2.9 - 10.5) (175/0.3)1/5 (/i/0.7)2/5 MeV. Remarkably, this is very sim- 

ilar to the range (mg = (6.6 — 15.4)r/^/^MeV) required in order for the self- 
interaction (r//4)(S'^S')^ to account for self- interacting dark matter when t] is 
not much smaller than 1. The corresponding coupling is As ~ (2.7 x 10"^*' — 
3.6 X 10^^) (175/0.3)2/5 (/^/0.7)^/^ implying that such scalars are very weakly 
coupled to the Standard Model sector. 
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1 Introduction 



It has become apparent that conventional colhsionless cold dark matter (CCDM) may 
have problems accounting for the observed structure of galaxies. N-body simulations 
with CCDM indicate that galaxies should have singular halos § with large numbers 
of subhalos p, Observationally, the density profile of galaxies in the inner few 
kiloparsecs appears to be much shallower than predicted by numerical simulations (the 
central density of dark matter halos being 50 times smaller than the CCDM prediction 
for dwarf galaxies and roughly independent of halo mass whilst the number of 

dwarf galaxies in the Local Group is an order of magnitude fewer than predicted 0, 
In addition, the CCDM predictions for the TuUy-Fisher relation |^ and the stability 
of galactic bars in high surface brightness spiral galaxies are not in agreement with 
what is observed, indicating lower density galaxy cores than predicted by CCDM. 
Although there is at present considerable uncertainty regarding the interpretation 



of observations and simulations 0, 0, it has nevertheless been argued that all the 
discrepencies between observations and simulations may be understood as indicating 
that dark matter halos in CCDM simulations are more centrally concentrated than 



observed 1 1 



In order to overcome the possible deficiencies of CCDM halos, one suggestion has 
been that the cold dark matter particles have a non-dissipative self-interaction P^ , 
and it has been shown that such cold, non-dissipative self-interacting dark matter 



(SIDM) can be effective in alleviating the various problems of CCDM Scattering 
of dark matter particles stops gravitational accretion at the centre of the halo and 



so allows a smooth core to form. Simulations with SIDM [0, are able to simul- 
taneously account for the observed density profiles of galactic halo and the number 
of subhalos [O. In the future SIDM may be strongly constrained by gravitational 



lensing observations of the shape of cluster halos [|T^, |T^, 17] and by the formation of 
massive black holes at the centres of galaxies, which is enhanced by self-interactions 
of dark matter particles |TB[ . 



In order to be able to account for the observed properties of dark matter halos. 
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the requirement on the mass M and self-interaction scattering cross-section a of the 
SIDM particles is that Pl 



rs = ^ = (2.05 X 10^ GeV - 2.57 x 10^) GeV-=^ . (1) 

The upper bound corresponds to the limit at which galaxy halos in massive clusters 
are destroyed by interacting with hot particles in the cluster halo (evaporation) [|13| 
whilst the lower bound corresponds to the limit where the SIDM particle would not 
interact within a typical galactic halo during a Hubble time [|12|, |l^ . 

The canonically simplest dark matter particle is arguably a gauge singlet scalar S. 
The possibility that gauge singlet scalars, interacting with the Standard Model sector 
via a coupling to the Higgs doublet of the form S'^SWH, could naturally constitute 
dark matter has been pointed out by a number of authors in the past ^ as well as 
more recently These calculations consider the case of massive (> 1 GeV) scalars 



which freeze out of thermal equilibrium when non- relativist ic p2|]. However, the range 



of 5* masses considered is too large to account for SIDM with perturbative S couplings. 

It has recently been noted that gauge singlet scalars have a natural self-interaction 
via an S'^-type coupling and so in principle could account for SIDM ^ An 
estimate of the upper bound on the coupling of S scalars to the Higgs doublets for 
5" mass of the order of 10-100 MeV (which is of the greatest interest in the case of 



perturbative S self-interactions) was derived in [|2J] by requiring that S scalars do not 
come into thermal equilibrium and so overpopulate the Universe. 

In this letter we consider thermal generation of SIDM S scalars which do not 
achieve equilibrium. We will show that such non-equlibrium thermal generation can 
naturally account for a dark matter density of 5* scalars with the right properties to 
account for SIDM. 

The letter is organized as follows. In Section 2 we discuss the perturbative upper 
limit on the S scalar mass. In Section 3 we consider the thermal generation of a relic 
density of S scalars. In Section 4 we consider case of zero bare S mass and the resulting 
consistency of the relic density, S mass and Spergel-Steinhardt SIDM cross-section for 
natural values of the S self-coupling. In Section 5 we present our conclusions. 
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2 Limit on ms for Perturbative SIDM 



We first consider the perturbative upper limit on the S mass if it is to play the role of 
SIDM. We will consider the case of complex gauge singlet scalars for consistency with 
the cross-sections and discussion given in [H, which we will use here. We expect that 
the results for real scalars will be very similar. The model is described by 



The total centre-of-mass 5* scattering cross-section is the sum of SS'' 

ss ss, 

= c^sst^sst + c^ss^ss 



1287rm| ' 



Therefore 



which implies that 



rs 



ms 1287rm| ' 



(2) 

SS'' and 
(3) 

(4) 



ms = 35.8aJ/^ 



1/3 /2.05 X 10^ GeV-^y^^ 



MeV 



rs 



(5) 



where = rf /4:Ti. (Similar expressions have been obtained in |]2^ , ^.) Thus if we 
require that a,, 1 in order to have a perturbative theory then the condition Eq. (|l|) 
requires that ms = 0^^(15.4 — 35.8) MeV ^ 30 MeV. (We refer to this as the Spergel- 
Steinhardt mass range.) We note that this puts a severe bound on the coupling A5, 
since the S scalar gains a mass from the Higgs expectation value, 



mc 



m^ + 



(6) 



where v = 250 GeV. Thus the requirement that ms ^ 30 MeV imposes an upper 
bound on A5 



Xs < 2.9 X 10" 



ms y 
30 MeV/ 



(7) 



The Spergel-Steinhardt mass range assumes a perturbative 5* self-coupling. A non- 
perturbative self-coupling may be possible, but it would be difficult to calculate the 
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properties of such a model, so we must restrict ourselves to the perturbative case. In 
addition, the only known scalar self-coupling, that of the Standard Model Higgs dou- 
blet, \h{H^HY, is given by = m\/Av'^ = 0.053(m/i/115 GeV)^. (The experimental 
lower bound on the Higgs mass is rrih > 113 GeV I^G], whilst an upper bound for the 
pure Standard Model is obtained from radiative corrections to electroweak observables. 



rrih < 165 GeV ||2^. The upper bound in extensions of the Standard Model can be 1 



TeV or larger p8|.) This is typically perturbative but not very much smaller than 1, 
suggesting that a natural value for the S scalar self-couplings is around 0.1. 



3 Thermal Generation of S Scalars 

There are two processes which can produce a density of S scalars: 2^2 annihilation 
processes and decay of a thermal equilibrium density of Higgs scalars to SS"^ pairs, 
h° S'S'T. We first consider 2^2 annihilations. The relic density from scatter- 
ing processes in a radiation dominated Universe is found by solving the Boltzmann 



equation |2^, |2£ 



- -\j - Jo) ; ^ 



(8) 



dT K 

where / = ns/T^, fo = n-o/T^ and g(T) = Qb + Tqf/^, where qb and gp denote the 
number of relativistic bosonic and fermionic degrees of freedom respectively, is the 
number density of S scalars and rio is the thermal equilibrium S number density; for 
relativistic 5* scalars 

= (^) . (9) 

We consider the case where the S scalar density is very small compared with the 
equilibrium density and solve the Boltzmann equation with / = on the right-hand 
side, 

dT ~ K ^° ■ ^ ' 

We take the Universe to be initially at a high temperature, T ^ mw, and calculate 
the resulting relic density of 5* scalars as the Universe cools. The annihilation cross- 
sections of relativistic SS'^ pairs to t quarks, W and Z bosons and the h° Higgs 



scalars (lighter quarks and leptons do not contribute significantly due to their very 
small Yukawa couplings) are estimated by using the centre-of-mass annihilation cross- 
sections calculated for 5* scalars with typical energy Et ~ T . We will see that the 
core results of the paper are not very sensitive to uncertainties in the calculation of 
the annihilation cross-section and thermal relic S density. The relativistic annihilation 



cross-sections ai may be obtained from the non-relativistic < aannVrei > given in 
via the relation CTj = (1/2) < (TannVrei > (^5 — ^ Et) (where i = t,W,Z,h" denotes the 
Standard Model particle in question), which may be confirmed by directly calculating 
the cross-sections. For T < rrii the contribution of o"j to the total cross-section is zero, 
which models Boltzmann suppression. Then < aiVrei >= 2o"i, where we take Vrei = 2 
for relativistic annihilations [^. For relativistic 5* scalars, fo = 1.2/ -k"^ is a constant 
so Eq. (0) can be integrated as 



h = -VI dEr ^ , (11) 



where Et = T, Etj = rrii and the initial thermal energy Et^ — > 00. We will take 
K oc g{TY^'^ to be constant with g{T) = g(Ti), where Tj = rrii, since most of the 
integral comes from Et close to rrii. The total contribution to / is then 

In this we have taken g(Ti) = 106.75, corresponding to the Standard Model degress 
of freedom and At = 0.7 (corresponding to rrit = 175 GeV). In addition, in order to 
obtain an analytical result we have expanded the Higgs propagators in cTj {i = W, Z, t) 
assuming that 4i?|. is large compared with m\, which is generally satisfied if m\ is 
small compared with ^m\r. (We refer to this as the small Higgs mass limit.) 

The S number density from the decay of thermal equilibrium h° scalars at temper- 
atures less than the electroweak phase transition (where Tew ^ 1.5mh [0) is given 

by 

+ 3Hns =< Th" > riho , (13) 
where H is the expansion rate and the thermal equilibrium density of h°, rih" eq, is 
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given by 

1 E {E'' - ml^^^^ 



and the decay rate for h° scalars with energy E is 

r.o = ^ . (15) 

Thus the thermal average of the decay rate is 

Therefore in terms of /, the S density from h° decays is given by 

df_ ^ <Tho>fhoe, ^ viniH/T) Xlv'e-^'^/^ 
dT KT^ ~ KT^ 327r3 



:i7) 



ri{a) is a slowly varying function of a, with 77(0) = 1.64, 77(1) = 1.87 and 77(5) = 3.00. 
Since most of the contribution to / comes from rrih/T ~ 1, we take ri{mh/T) to be 
equal to ^^(l), in order to obtain an analytical expression. Therefore the density of S 
scalars from h° decay, fdec, is given by. 

In this we have assumed that v is given by its T = value, v = 250 GeV. Since most 
of the contribution to fdec comes from T ^ rrih < Tew, this should be a reasonable 
approximation. We see that the S density from h° decays is generally much larger 
than that from 2 2 annihilation processes in the small Higgs mass limit. For larger 



values of the Higgs mass, it is possible that s-channel pole annihilations may result 



in 2 ^ 2 processes dominating the h° decays |31|, in which case fdec is a lower bound 
on the number of 5* scalars produced thermally. 

The resulting density of S plus scalars is then the sum of scattering and decay 
contributions 

where pc = 7.5 x 10~^'^h'^ GeV^ is the critical density, = 2.4 x 10~^ eV is the present 
photon temperature, Tj ^ rrii and g{Ty) = 2. In this we have used the fact that the 



S number density to entropy is conserved once the scattering and decay processes are 
Boltzmann suppressed, such that g(T)ns/T^ is constant. Therefore with / fdec and 
g(Ti) = 106.75 for Tj f» rrih, the thermal rehc S density Qs is related to \s by 



A. 



Thus for Higgs masses in the range 115 GeV to 1 TeV and with expansion rate h 0.7, 
the upper bound on A5 from requiring that Vis ^ 0.3 is in the range (1.4 x 10^^° — 
3.6 X 1^-^)7]-^'^. This is in broad agreement with the upper bound estimated in 
based on the weaker condition that the S scalars do not come into thermal equilibrium. 
More importantly, we see that it is possible to generate a thermal relic density with 
Vis ~ 0-3 and 10 MeV (typical of SIDM scalars) purely within the minimal 

gauge singlet scalar extension of the Standard Model. 



4 Naturally Consistent Thermal Relic SIDM for 
Zero Bare Mass 

The value of A5 from requiring that Vis ~ 0.3 is satisfied is not very much smaller 
than the upper limit Eq. (0) from the requirement that the Higgs expectation value 
contribution to the S mass is compatible with perturbative SIDM S scalars. This 
suggests that it is quite likely that all the 5* mass might come from its interaction with 
the Higgs scalar when its relic density is sufficient to account for dark matter. If we 
assume that all the 5* mass is due to the Higgs expectation value, then we find that 
the S mass is fixed by the thermal relic density 

„,.2.9f^k)"7±y'V^)'%,eV. (21, 
* V0.3/ V115GeV/ ^ ' 

We refer to this as the thermal relic S mass. Comparing with the Spergel-Steinhardt 
range for SIDM, 

ms = (6.6 - 15.4)77^/3 MeV , (22) 

we see that the thermal relic mass for S scalars is within the range required to account 
for SIDM when the self-coupling constant rj is equal to about 0.1, a natural value 



which is consistent with the Higgs doublet self-couphng in the Standard Model. The 
thermal relic mass is not strongly dependent upon cosmological parameters, nor is it 
strongly dependent upon the Higgs mass. In particular, it is relatively insensitive to 
uncertainties in the calculation of /, since a change in / by a factor S produces a 
change in Qs by the same factor, and so a change in the thermal rehc mass by 6^^^. 
The coupling corresponding to the thermal relic mass is 

X,^2,.10-'i^f(jLr(^f. (23) 
V0.3; \0.7 J Vll5GeV; ^ ^ 

This suggests a scenario for dark matter in which stable gauge singlet scalars couple 
very weakly to the Standard Model sector but self-couple with a relatively strong 
coupling of about 0.1, of the order expected from the example of the Standard Model 
Higgs self-couphng. 



5 Conclusions 

We have considered the thermal generation of a relic density of self- interacting dark 
matter (SIDM) gauge singlet scalars. The dominant process for small Higgs mass 
is the decay of thermal equhbrium Higgs scalars to gauge singlet scalar pairs. For 
SIDM scalars with perturbative self-interactions, the mass must be no greater than 
around 30 MeV. For such light scalars, the requirement of an acceptable relic den- 
sity of S scalars requires that the S coupling to the Standard Model Higgs satisfies 
Xs ^ 10~^^~^^\ This limit comes from the requirement that S scalars are not thermally 
overproduced, which is a stronger condition than requring that they do not come into 
thermal equilibrium. 

In the case where S scalars account for dark matter and where the S mass is 
entirely due to the Higgs expectation value, wc find that the S mass is fixed by the 
thermal relic dark matter density to be between about 2.9 MeV and 10.5 MeV for 
Higgs masses ranging from 115 GeV to 1 TeV. (The upper limit on the S mass may be 
smaller if 2 2 annihilations dominate h° decays for large Higgs mass.) This is very 
similar to the range of masses ( (6.6 — 15.4)?7^/"^MeV) required by self-interacting dark 
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matter with self-coupling t] of the order of the natural value (based on comparison 
with the Standard Model Higgs doublet self-coupling) of around 0.1. This result is 
not strongly sensitive to uncertainties either in the cosmological parameters or in the 
calculation of the thermal relic 5* density. We find this coincidence remarkable and a 
possible hint that light gauge singlet scalars with very weak coupling to the Standard 
Model sector may play an important role in cosmology and particle physics. 

The Author would like to thank the CERN theory division for its hospitality and 
PPARC for a Travel Fund award. 
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